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Abstract 



.^ I The Hubbard model with an additional bond-charge interaction X is solved 

V^ I exactly in one dimension for the case t = X where t is the hopping ampli- 

^— ^ ■ tude. In this case the number of doubly occupied sites is conserved. In the 

sector with no double occupations the model reduces to the U = oo Hubbard 
'^ , model. In arbitrary dimensions the qualitative form of the phase diagram is 

obtained. It is shown that for moderate Hubbard interactions U the model 
C^ I has superconducting ground states. 
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I. INTRODUCTION 

The one- dimensional Hubbard model is since long the prototype of an exactly solvable 
model for correlated electrons |jl|,0| . However, it is difficult to include additional interactions 
so that the resulting model is still integrable. In this paper we present a generalized Hubbard 
model which apart from the Coulomb interaction U also contains a bond-charge interaction 
term X. This model has been studied extensively by Hirsch ^j^ who argued that it might 
be relevant for the description of high-T^ superconductors. For certain values of X and large 
densities of electrons (small doping) the bond-charge interaction can lead to an attractive 
effective interaction between the holes and the formation of Cooper-pairs. Although this 
picture was confirmed by a BCS-type mean-field theory it is desirable to find exact results 
confirming this behaviour. A first step in this direction was made in |^ where a simplified 
version of Hirsch's Hamiltonian has been studied in one dimension by Bethe-Ansatz methods. 
Indeed, it has been found that this model has a strong tendency towards superconductivity 
(see also 0). 

On the other hand, there has recently been a lot of activity in the investigation of the 
so-called r^-pairing mechanism of superconductivity. This idea has been introduced by Yang 
for the Hubbard model. It allows for the construction of states exhibiting off-diagonal- 
long-range-order (ODLRO). As shown in [BUI] ODLRO implies also the Meissner effect 



and flux quantisation and can thus be regarded as deflnition of superconductivity. In the 
case of the Hubbard model it was found that the ground state is not of the r^-pairing type 
0. The flrst model with a truly superconducting ground state of the //-pairing type is the 
supersymmetric Hubbard model introduced by EBler, Korepin and Schoutens (EKS model) 



11,12]. The EKS model is a Hubbard model which contains additional nearest-neighbour 
interactions. It has some attractive features. Apart from being exactly solvable by Bethe 
Ansatz in one dimension it is possible to determine the T = phase diagram in arbitrary 
dimensions. One always flnds a superconducting ground state ii U < Uc where the critical 
value Uc is positive corresponding to a repulsive on-site Coulomb interaction. 



In [jT3[ it has been shown that an r^-pairing ground state is not an exotic phenomenon 
restricted to one special model but may be found for a large class of Hamiltonians. As a 
special case the Hubbard model with bond-charge interaction X is in this class, provided 
that X = t where t is the hopping matrix element. It is this model which will be studied in 
this paperQ 

In Sec. II the symmetries of the Hamiltonian are investigated. It turns out that the 
model under investigation is a generalization of the U = oo Hubbard model to a model with 
a conserved (but in general non-zero) number of doubly occupied sites. Since the U = oo 
Hubbard model has attracted a lot of interest in recent years J]T1|-^ it seems worthwhile to 
study any generalization of it. 

In Sec. HI we present the exact solution in one dimension. It will be shown that the 
Hamiltonian can be mapped in certain subspaces onto a spinless free fermion model with 
twisted boundary conditions where the twist depends on the subspace considered. Using 



^A brief account of some of the results presented here can be found in jl3[|. 



this fact it is easy to obtained the spectrum and the zero temperature phase diagram of 
the Hamiltonian. It is found that for U not too large the model contains superconducting 
ground states of the r^-pairing type. For certain densities this is true even for moderately 
repulsive Coulomb interactions. 

In Sec. IV we construct the qualitative form of the phase diagram in arbitrary dimensions. 
Basically it has the same form as the phase diagram in one dimension although the exact 
location of all phase boundaries can not be determined. 

Sec. V contains a summary of the results and a discussion of the stability of the super- 
conducting ground states. 

In the Appendix details of the exact solution in one dimension are given. 

II. THE MODEL AND ITS SYMMETRIES 

The Hamiltonian of the Hubbard model with an additional bond-charge interaction X (or 
correlated hopping) on ad dimensional lattice with L sites and periodic boundary conditions 
is given by 

L 

n{X, U) = -tY, H (c]aQ<x + clcj^) +UY, nj^nji 
m ^=T,i i=i 

+^mi [c]aCl^ + cLcja) {Uj^-a + ni^-a) ■ (2.1) 

Here Cja, Cj„ are the usual Fermi operators and rij^ = Cj^Cja is the corresponding number 
operator, (jl) denotes nearest-neighbour sites on the d-dimensional lattice. In the following 
we will be interested in the special case 'H{U) = 7i{X = t,U). In this case it is convenient to 
rewrite the Hamiltonian in terms of Hubbard operators X"* = \a){b\ (a, 6 = 0, ±1, 2) where 
denotes an empty site, ±1 a site occupied by a single electron with spin | or |, and 2 a 
doubly occupied site. The fact that for each site j the four states |0), | — 1), |1), and |2) 
form a basis of the local Hilbert space leads to the local constraint X^'^ + Xj"^ + J2o- ^'j'^ = 1- 
The Hubbard operators obey the following graded commutation rules 

[Xf, X[% = (xfS,, ± Xf5ad) 5,1 , (2.2) 

where the anticommutator has to be taken only if both operators are fermionic, i.e. change 
the particle number by one (e.g. XJ^ or -V|°"). The standard fermion operators Cja can be 
expressed through the Hubbard operators 

4 = Xf + aX|'-^ c,^ = X]" + aX-"\ n,, = XY^Xf (2.3) 

and vice versa, e.g. 

XJ-^ = (1 - n,,_>,., Xf = (l-n,J(l-n,T). (2.4) 

In terms of the Hubbard operators the Hamiltonian ( |2.1| ) takes the following form: 



n{u) =no + nu 

L 

= -tJ2T. {xfxj"" + xf °xO" + xfx^'' + xfx^^") + UY, xf . (2.5) 

The Hamiltonian 'H{U) has a lot of symmetries, since only two types of processes are 
allowed: (i) a particle from a singly occupied site may hop to an unoccupied neighbour site, 
or {ii) a particle with spin a from a doubly occupied site may hop to a singly occupied 
(particle with spin —a) neighbour site. Therefore it is clear that 'H{U) conserves not only 
the total number A^ of electrons but also the number A'^}'^ of single electrons with spin a 
and the number N2 = SJLi ^t'^iJ. ^^ doubly occupied sites, i.e. ['Ho,'Hu] = 0. In addition, 
Tio has two SU(2)-symmetries: apart from the SU(2) spin symmetry, i.e. Tio commutes with 
5+ = Ej^=iS79a, S- = (5+)^ S' = l{Ni^^ - nI^^), Ho also commutes with the following 
so-called //-operators: 

^ = Ec,TC,i, ^^ = EsV]^, v' = l{N-L) (2.6) 



r(T) , at(1)) 
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where A^ = A'^i + 2A''2 (with Ni = A^} + Ni ) is the total number of particles. In terms of 
the Hubbard operators the r^-operators are given by r^ = — J2j=i X^'^ and r]^ = — Z^JLi X'^'^. 
Tj^ creates a double occupation with momentum zero as can be seen from its form in mo- 
mentum space ?7^ = J2k CkiC^ki where c[^ = -^ J2k e~^-'''Cj^. Thus t]'^ is just the conventional 
s-wave pairing operator. Note that in the case of the Hubbard model the corresponding 
?7-symmetry is generated by 77,^ = Y,^=i{~^)'' Cj^Cj^ corresponding to pairs with momentum 

The Hamiltonian 7i{U = 0) without the Coulomb interaction is also invariant under a 
particle-hole transformation: 

Uc,M^' = cl, Uc],U-^ = c,., (2.7) 

UUdA-^ = Ho , (2.8) 

whereas the Coulomb interaction Tiu = U I]j=i i^ji^ji transforms as 

UHuU-^ =nu + U{L-N). (2.9) 

Due to this particle-hole symmetry we can restrict our investigation to the case N < L. 

In the following especially the conservation of N2 will be important. It allows to diago- 
nalize the Hamiltonian in subspaces with fixed A"2. For A^2 = the Hamiltonian Tio reduces 
to the ?7 = cx) Hubbard model ITl-El 



Hoc = -t E E {xfx^'' + xf °x°'^ 



= -^ ^E E (4c,. + 4q.) V (2.10) 

iji) '^=Ta 

where V is the projector onto the subspace with no doubly occupied sites. Thus the Hamilto- 
nian ( |2.5| ) is a generalization of the the U = 00 Hubbard model to a model with a conserved 



(but in general nonzero) number of doubly occupied sites. This is similar to the EKS-model 
||TT| , [T2[| which is a generalization of the supersymmetric t-J model [^ to a model with a 
conserved number of doubly occupied sites. 

III. EXACT SOLUTION IN ONE DIMENSION 

In the following the Hamiltonian TYq (see ( ^.51 ) ) is diagonalized in one dimension with 
periodic boundary conditions^]. 

In the subspace H^^^]^^ with Ni singly occupied sites and N2 doubly occupied sites, i.e. 
A^ = A^i + 2N2 particles, we define the following basis vectors: 

|f,a,6) = n X:f n Kf\^) (1 < a^i < ■ ■ ■ < x^, < L). (3.1) 

a=l /3=1 

|x, CT, 6) is a state in which the sites Xa are occupied by a single electron with spin Ua- 
{Vii ■ ■ ■ Vl-Ni} = {1; • • • ) -^} \ {^1; • • • ; X jy^^} are the sites occupied by a boson ba where 



ba = 2 for a doubly occupied site and &« = for an empty site. Note that due to ( |2.2| ) all 
Xy'^' commute mutually and with all the X^^'^. 

Next we define the operators C and C generating cyclic permutations of the spins and 
bosons, respectively: 

C{ai,. . . ,(TnA = W2,- ■ •jO^ATijO'i}, (3.2) 

C {61, . . . , ^Ar-AT-^} = {6Ar_Ar^_i, 61, &2, ■ • • , bN-Ni~l}- (3.3) 

Using these operators we can define subspaces of His[i,N2- For every spin configuration 
{(Ti, . . . , (Jni) there exists a minimal integer K >1 such that 

C^{ai,...,o-7vJ = {o-i,...,cr7vJ. (3.4) 

Clearly for a fully polarized ferromagnetic state we have K = 1 and for a Neel state K = 2. 
In the same way we define the integer K' for the distribution of bosons. 

All the states \x, a,b) with a and b characterized by the same integers K and K' span 
a subspace Hni,N2{K,K') of Hj^^^n^. It is now important to notice that the subspaces 
Hni,N2{K, K') are invariant under the action of the Hamiltonian Ti. This fact is well-known 
p!^ , p!6[] for the U = oo Hubbard model, i.e. N2 = 0. Since the local Hamiltonian /ijj+i only 



permutes bosons with fermions a spin cTq+i will under the action of H always stay "to the 
right" of a spin a^ except for the case where it moves "over the boundary" {L ^ L+1 = 1). 
We can now restrict ourselves to to the diagonalization of the Hamiltonian in the sub- 
spaces H]\f-^^^]s[^{K,K'). In the appendix it is shown that in each of these subspaces TC is 
equivalent to a free fermion Hamiltonian Ti with twisted boundary conditions: 



^A solution of an equivalent model with open boundary conditions has been given in [^ (see also 

1). 



L-1 

H = - ^ (4a,+i + 4+ia,) - {e'^^'^a{a^ + e-'^^'^ajai) (3.5) 



which now acts on the "stripped states" [0] 



AT 



}=n4jO). (3.6) 



lyy — 


L 




2n 


k = 


k' 


k' = 


K' 



The aj are (spinless) fermion operators and the allowed values A0 for the twist in the 
boundary conditions is different for each subspace: 

■ ' ^'~^ (3.7) 

(z/ = 0,l,...,ir-l), (3.8) 

(i.' = 0,l,...,ir'-l). (3.9) 

After the canonical transformation a] — *> e^^^'^a^^ aj — > e~^^^'^aj the Hamiltonian (|3.5|) is 
seen to be equivalent to the following translational-invariant free fermion Hamiltonian 

L 

7i:,g = - ^ (e^^*a]aj+i + e-^^'^aj+ia^) (3. 10) 

from which the eigenvalues of Ti can be obtained easily by Fourier transformation: 

E = -2^cos(g, + A<^)n,^ (3.11) 

V 

with I^jy^?,. = -^1 and Uq^ = 1 (0) if the mode q^ is occupied (not occupied). The wavenum- 
bers q^ can take the values q^ = ^u (z/ = —^ + 1, . . . , ^). 

In the ground state the g^, are as symmetric as possible around z/ = 0. For A^^i even we 
have Ug^ = 1 for z/ = —^ + 1,...,^ and thus 

cos f f + A0) 

^ = -2 .^, ,M ^ sin(iVi7r/L) . (3.12) 

sm (vr/L) 

For A''i odd we have to choose u = —^^^ + 1, . . . , ^^^^ yielding 

£._2£2!iM,i„,Ar.,/i). (3.13) 

sm (tt/Lj 



This shows that in the ground state one has 



—n/L {Ni even), 
(A^i odd), 



(3.14) 



which gives the ground state energy 



E, . _,M^, (3.15, 

sm (vr/L) 

The ground state energy only depends on the number A^i of singly occupied sites and is 
independent of the number A^2 of doubly occupied sites. Therefore the ground state energy 
of (|3D is Eo{N,, N2, U)=Eo + UN2. 

The mapping of the original model onto spinless free fermions ( ^.10| ) implies that the 
eigenvalues of (|2.5|) are highly degenerate in general. 



IV. PHASE DIAGRAM AT T=0 



In this section we determine the phase diagram of the Hamiltonian ( p.5| ) in the U-D- 
plane where D = N/L is the particle density. Firstly, we consider the one-dimensional case 
which can be treated exactly using the results of the previous section. Then we show that 
the phase diagram in higher dimensions has qualitatively the same form as in one dimension. 
We find the same phases although we cannot determine all phase boundaries exactly. Finally 
we calculate correlation functions for generalized ?7-pairing states which appear as ground 
states in certain parameter regions. 

A. Phase diagram in one dimension 

In order to determine the phase diagram in the [/-D-plane for the one-dimensional case 
we have to minimize the ground state energy^ Eo{Di, D2, U)/L = — - sin{Di7T) + UD2 for a 
fixed particle density D = Di + 2D 2- A simple calculation yields 

ro (t/<-4), 

Di = I i arccos(-t//4) (-4 < f/ < [7^), (4.1) 

[d iu>u,), 

where Uc{D) = -4cos(7r£)). 

Due to the particle-hole symmetry it is sufficient to discuss the phase diagram for D < 1. 
We find four different phases which will be discussed separately in the following. 
Regime I : U < —4 

From ( [4.1D we see the ground state contains only doubly-occupied sites and no single elec- 
trons. In this case the ground state energy is simply Eq{N, 0,U)/L = UN/2. In the absence 
of single electrons the double occupations are static and all states with the same number 
A^2 = N/2 of doubly-occupied sites have the same energy. Among these states are the 
generalized ?7-pairing states 

\^p) = ivl^f^'lO), (4.2) 

4 = E e'^'^sV], (4.3) 



^Froni now on we work in the thermodynamic limit L, Na — > 00 with Da = N^/L fixed (a = 1, 2). 



with momentum P = ^z/ (z/ = 0, 1, . . . , L — 1). These states are ground states for all chain 
lengths L, not only in the thermodynamic limit. All of these states show ODLRO, i.e.0 

As a consequence these ground states are superconducting since it has been shown that 



ODLRO also implies the Meissner effect and flux quantisation [p|-[TO| . 

Regime II : -A < U < Uc{D) = -Acos{ttD) 

For —4 < U < Uc{D) the ground state has both a finite density of single electrons and 

of double occupations. Again it is highly degenerate. Here the ground state energy is 

Eq{Di,D2,U)/L = — 2^Vl6 — f/^ + UD2. Among these ground states are the ?7-states 

(?7o) ^\U = oo) where \U = oo) stands for an arbitrary ground state of the U = oo Hubbard 

model at particle density Di = - arccos(— f//4) and N2 is then obtained from D2 = N2/L = 

[D — Di)/2. Again these states have ODLRO (see the discussion of correlation functions 

below) and are thus superconducting. It is interesting that this superconducting phase 

extends into the region of positive U, i.e. Coulomb repulsion. 

Regime HI : U > U^iD) = -A cos{7rD) and D < 1 

In regime III we have no doubly occupied sites in the ground state and every ground state of 

the U = 00 Hubbard model is therefore also a ground state of (4). Regime HI' with D > 1 

is obtained by using the particle- hole symmetry ( |2.7| )-( |2l9| ). 

Regime IV : U > Uc{D) = -A cos{7rD) and D = 1 

This case is much like regime HI, but now the ground states \U = 00) are insulating. The 

point (-D = 1, f/ = 4) corresponds to a metal-insulator transition. 

The complete phase diagram for the one-dimensional case is shown in Fig. 1 (see also [^,^). 

B. Phase diagram in arbitrary dimensions 



In higher dimensions the phase diagram can be constructed along the lines of ||T2[. In 
13| it has been shown that \ipo) = (?7q) |0) will be a ground state for U < —2Z where Z 



is the number of nearest neighbour sites in the d- dimensional lattice. In order to construct 
the full phase diagram we need the following three properties: 

1) r^-symmetry, i.e. [^0,^70] = 0, 

2) conservation of the number N2 of doubly occupied sites, 

3) for f/ = the ground state energy does not depend on the number N2 of doubly occupied 

sites. 

The first two properties have already been demonstrated in Sec. II, property 3) can be proven 
by generalizing the argumentation of ||2^ where the analogous property for the EKS-model 
has been derived. 



*See also Eq. (4.5) below. 



Using l)-3) the qualitative form of the phase diagram can be estabhshed in complete 
analogy with |12| by first considering the grand canonical ensemble and then translating the 
results into the canonical ensemble. One finds a phase diagram which looks very similar to 
that of the one- dimensional case (Fig. p. The same phases appear, only the location of the 
phase boundaries changes. Except for the boundary between regimes I and II (see |T^) we 
have not been able to determine them exactly. An interesting open question is thus whether 
the superconducting regime II extends into the positive-?/ region in all dimensions. 

C. Correlation functions 

Correlation functions with respect to the //-pairing states can be calculated in a straight- 
forward way 1^. Let |0) be a state with ?7p|0) = and r/^10) = ^{Ni — L)\(f)), e.g. |0) 



is a state with A^i singly and no doubly occupied sites. In our case we have |0) = |0) or 
10) = \U = oo). For the r^-pairing states \'ijjp{N2)) = {7]p) ^|0) we denote the expectation 
value of an arbitrary operator O by {O) = t^ (^^ iiI (^ )\ • We now can express any corre- 
lation function with respect to \iIjp{N2)) through correlation functions with respect to |0). 
With (C)^ = ^^g^ we find (for L^oo, D^ = N^/L and j ^ /): 

(4.Q.) = ^ \^_'j^^^" {clcir), - ^-^^''"-'-'^Y^l <^^^^^'-)^ ' (^-6) 

D2 ' 



+ [r^J ((i-^o(i-%))0, (4.7) 

{SfSn = {SfSn^ {a = x,y,z). (4.8) 

These results hold in all dimensions in the thermodynamic limit. Correlators for phase I 
(where |0) = |0)) are trivial. The asymptotics for the correlation functions in phase II 
(where |0) = \U = oo)) in one dimension can be obtained using the results of pB|-pT|. The 
result for the ?7-correlator {i]jr]i) proves the existence of ODLRO since the value of the limit 
in (pi) is found to be e*^('-^)(l - Di - £'2)^2 ^ 0. 



V. CONCLUSIONS 

In this paper a Hubbard model with an additional bond-charge interaction X has been 
investigated at the special point X = t where the number of doubly occupied sites is con- 
served. In one dimension the complete spectrum of the Hamiltonian could be determined 
exactly by mapping onto a system of spinless free fermions with twisted boundary conditions. 
In arbitrary dimensions the T = phase diagram in the D-U-p\ane has been obtained. One 
finds four phases, two of them containing ground states of the 77-pairing type. These states 
exhibit ODLRO and thus are superconducting. Therefore the model ( pTT|) provides a new 

9 



example of a purely electronic model with a superconducting phase. The Cooper pairs have 
zero-size which might be regarded as an approximation to the small coherence length found 
experimentally in high-T^ superconductors^. 

Let us finally discuss the question of stability of the superconducting phases. As in the 
EKS model all interaction constants (apart from the Coulomb interaction U) have to take 
certain values in order to allow for the //-pairing ground states. Especially the bond-charge 
interaction X has to be equal to the hopping matrix element t in order to guarantee the 
conservation of local pairs. One might wonder if the superconducting properties will survive 
if one allows for a decay of Cooper pairs (i.e. X ^ t). The results obtained in [^ suggest 
that this is indeed the case. 



For the EKS model other perturbations of coupling constants have been investigated ^5 
It has been found that in one dimension these perturbations destroy superconductivity, but in 
higher dimensions the supercondcuting phase is likely to be stable under such perturbations. 

A similar analysis is also desirable for the model presented here. It would be interesting to 
find a perturbation by allowing for (small) additional nearest-neighbour interactions which 
lifts the large ground state degeneracy in such a way that only (some) of the ry-pairing states 



remain as ground states. In p6|j37|] the ground state for the model ( |2.1|) has been obtained 
in the presence of an additional nearest-neighbour Coloumb interaction V at half-filling in 
arbitrary dimensions. In that case the ground state is unique (apart from a trivial twofold 
degeneracy) for U < 2ZV — Zmax(2t, V), but not superconducting. 

These questions are currently under investigation using perturbation theory and exact 
diagonalizations of small systems. Results will be presented in a future publication. 
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APPENDIX: 

In this appendix we extend the method used in [|l^] to map the U = oo Hubbard model 
onto a free fermion Hamiltonian with twisted boundary conditions to our case. 

First we investigate the action of the local Hamiltonian (we set t = 1 from now on) 

h _/,(+) ,/,(-) /,(-) _ //,(+) V 

^Sli = - E E xfx^i, (Ai) 

a=T,i 5=0,2 

on the states |x, a, b). We restrict ourselves to the investigation of /ijj+i since h)~j^-^ can be 
treated analogously. 



^For reviews on other models of local-pair superconductors, see e.g. [0;33] and references therein. 
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It is easy to see that hjj_^_^\x, <?, 6) = if either j G {xi , . . . , xni } or j + 1 ^ {xi, . . . , xn^}. 
Therefore it is sufficient to consider the case j + 1 = x-y E {xa}, j = yy G {Vis}- A 
straightforward calculation yields for 1 < j < L — 1 



,(+) 



with 



and for j = L 



with 



h'^/^,\x,a,b) = -\x',a,b) (A2) 



^«h;° \:ti] (A3) 



/ig|f,c?,6) = {-lf^\x~^',C{a^},C'm) (A4) 



„_rx„+i (a = l,...,A^i-l) , . 

In the subspace Hn-^^n^{K,K') we introduce the following states: 

\x, k, k') = Y. E e'^'e-''''^\x, C'{a4, C""{6^}) (A6) 

1=0 m=0 

where fc and k' are given by ( p.8|) and ( p.9|) , respectively, /ij^+i acts on these spaces in a 
simple manner: 

/,(+) 1^ L L/\ _ /-|2^',/^,A;') (j = i,...,L-i) ..„x 

in the case j + 1 G {xq,} and j G {y/3} and is zero otherwise. Thus Z^JLi ^jj+i ^^ts on the 
states |x, A;, A;') in the same way as — Yj^zI ctjCtj+i — e^^^'^a^ai acts on the "stripped states" 



( p.6|) where the a^, a^ are spinless fermion operators and A0 is given by 
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FIGURES 

FIG. 1. T = phase diagram in one dimension. The phase diagram in higher dimensions has 
quahtatively the same form. 
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